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Abstract 



We present density matrix renormalisation group calculations of the Pariser- 
Parr-Pople-Peierls model of linear polyenes within the adiabatic approxima- 
tion. We calculate the vertical and relaxed transition energies, and relaxed 
geometries for various excitations on long chains. The triplet (1 3 B^) and 
even-parity singlet (2 1 A^) states have a 2-soliton and 4-soliton form, respec- 
tively, both with large relaxation energies. The dipole-allowed (id-B") state 
forms an exciton-polaron and has a very small relaxation energy. The relaxed 
energy of the 2 1 A+ state lies below that of the l 1 ^" state. We observe an 
attraction between the soliton-antisoliton pairs in the 2 1 J 4+ state. The cal- 
culated excitation energies agree well with the observed values for polyene 
oligomers; the agreement with polyacetylene thin films is less good, and we 
comment on the possible sources of the discrepencies. The photoinduced ab- 
sorption is interpreted. The spin-spin correlation function shows that the 
unpaired spins coincide with the geometrical soliton positions. We study the 
roles of electron-electron interactions and electron-lattice coupling in deter- 
mining the excitation energies and soliton structures. The electronic interac- 
tions play the key role in determining the ground state dimerisation and the 
excited state transition energies. 

PACS numbers: 71.10.F, 71.20.R, 71.35 
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I. INTRODUCTION 



The inter-play of electron-electron interactions and electron-lattice coupling in linear 
polyenes results in a wealth of low-lying excitations. Electron-electron interactions induce 
spin density wave correlations in the ground state. The lowest lying excitations are triplets, 
which combine to form dipole-forbidden singlet ( 1 ^4jj~ ) excitations. Optical excitations are 
gapped, lie above the ^A^ state, and are essentially ionic in character, that is, there is 
charge transfer from one site to another. The lowest optically allowed (l 1 ^") state lies 
below the charge gap [Q, and is thus excitonic in character. For convenience, we show the 
group theoretic labelling of the states discussed in this paper in Table I. 

Electron-phonon interactions result in a dimerised semiconducting ground state. Within 
the adiabatic approximation, the non-linear excitations include charged-spinless and neutral- 
spin 1/2 solitons. Both electronic interactions and electron-lattice coupling lead to a gap in 
the optical spectrum. In contrast to the interacting limit, however, the 2 1 v4+ state always 
lies above the ^B~ state in the non-interacting electron-phonon model. 

The realisation that electronic interactions play a significant role in polyenes came via 
the experimental observation, by Hudson and Kohler in 1972, that the 2 1 A+ state lies 
below the l l B~ state. At the same time, by perfoming a double configuration interaction 
calculation on the Pariser-Parr-Pople model, Schulten and Karplus [[| demonstrated that 
the 2 1 A+ wavefunction has a strong triplet-triplet contribution, and has a lower energy than 
the l l B~ state. The triplet-triplet and correlated nature of the 2 1 A+ state has been further 
investigated by Tavan and Schulten 0] and other workers ||. In 1986, Hayden and Mele 
H performed a real space renormalisation group calculation on the Hubbard-Peierls model 
of up to sixteen sites and found that the 2 1 v4+ state was composed of 4-solitons. This 4- 
soliton nature has also been investigated by Su fjj, and Wen and Su [§]. Ovchinnikov et 
al. also high-lighted the role of electronic interactions, by suggesting that they are largely 
responsible for the optical gap ||. In contrast to the strong deviations from the ground 
state geometry predicted for the triplet and 2M+ state, Grabowski et al. [|ll| predicted that 
the l 1 ^" state is an exciton-polaron. 

The existence of the 2 1 A+ state below the l l B~ state in polyacetylene thin films has 
been suggested by a number of experiments. Third harmonic generation (THG) and two 
photon absorption by Halvorson and co-workers JTT] indicate that a 1 A+ state lies below 1.1 
eV; while the linear absorption, locating the ^B~ state, typically rises at 1.8 eV and peaks 



at 2.0 eV [ 121 . However, Fann and co-workers |T3[ performed THG, finding peaks at 0.6 eV 
and 0.89 eV, which they interpret as 1 A+ and 1 B~ states virtually coincident at 1.8 eV. The 
position of the 2M+ state is therefore not definitively established. We return to this point 
in section V when we discuss our own theoretical predictions. For a detailed review of the 
experimental and theoretical studies of conjugated polymers up to 1992, see |L4[ . 

Electron-electron interactions in 7r-conjugated systems, such as irons-polyacetylene, are 
conveniently modelled by the one-band Pariser-Parr-Pople model, which includes long range 
Coulomb interactions. This semi-empirical model has been extensively used to study the 
excited states of small conjugated molecules with a remarkable degree of success Jl5|. The 
Peierls model describes the electron-lattice coupling in the adiabatic limit. Thus, the Pariser- 
Parr-Pople-Peierls model is a realistic and accurate model of ^-conjugated systems, which 
captures their essential physics. In an earlier paper |L6| we performed accurate calculations 
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on this model using the infinite lattice algorithm of the density matrix renormalisation group 



(DMRG) method [17], [18]. The Hellmann-Feynman theorem was used to calculate the low- 
lying excited states and the lattice geometry associated with them. We showed that the 
1 3 B+ and 2 1 /1+ states are modelled by 2 and 4 soliton fits, respectively, and that the l 1 -^" 
state is an exciton-polaron. In this paper we develop that work. In particular, our objectives 
are: 

1. To further demonstrate that the DMRG calculations are reliable by, (i) making compar- 
isons to the exact non-interacting limit, and (ii) comparing the infinite lattice method 
to the finite lattice method. 

2. Use a realistic model of polyenes to understand the roles of electron-electron interac- 
tions and electron-lattice coupling in determining the dimerisation of the ground state 
and the transition energies of the excited states. In agreement with the earlier work of 
Horsch [Tj|, and Konig and Stollhoff [20|, we find that the electronic interactions play 



the key role in driving the ground state dimerisation. Electronic interactions are also 
dominant in determining the solitonic structures and transition energies of the excited 
states. 

3. To make more detailed comparisons to other experimental probes, in particular photo- 
induced absorption. The agreement with a wide range of experiments confirms the 
validity of the model, our calculational method, and our predictions on the soliton 
structures and their interactions. 

4. To further investigate both the geometry and electronic properties of solitons. 



This paper also serves as a correction to |Tf| . In that paper we used the dimerised ground 



state geometry in the Coulomb interactions to calculate the energy of all the states. Thus, 
the Coulomb interactions (unlike the one-electron transfer integrals) were not updated in 
the Hellmann-Feynman minimisation procedure for the relaxed states. We find that using 
the correct geometry in the Coulomb interactions affects the excitation energies by ca. 0.1 
eV. The geometry of the triplet excited state is modified, so that now there is no soliton- 
antisoliton confinement in the triplet state. However, attractive soliton interactions remain 
in the 2 1 A+ state. 

The plan of this paper is as follows. In the next section we introduce the Pariser-Parr- 
Pople-Peierls model. To establish the consequences of the inter-play of electron-electron 
interactions and electron-lattice coupling, we consider these two limits separately in sections 
III and IV. The non-interacting limit also allows us to compare the infinite and finite DMRG 
algorithms to an exact calculation. In section V we solve the full model, and discuss the 
vertical and relaxed energies of the key excited states. As well as linear absorption and 
non-linear optical spectroscopies, photo-induced absorption is a useful tool in determining 
the positions of excited states. We discuss the experimental situation and our theoretical 
interpretation in section VI. In section VII we consider the solitonic structures. By mak- 
ing comparisons between the geometrical soliton structures and the spin-spin correlation 
functions, we show how they are closely related. We conclude and discuss in section VIII. 

As well as the work already mentioned, earlier work on the solitonic structure of the 



low-lying excitations include, a mean-field study of the Heisenberg-Peierls model [21] and an 
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exact diagonalisation of a 12 site extended Hubbard- Peierls model [22]]. The DMRG method 



has recently been used by Yaron et al. |23|] and Fano et al. [24| to solve the Pariser-Parr-Pople 
model for linear and cyclic polyenes, respectively. 



II. THE PARISER-PARR-POPLE-PEIERLS MODEL 

The Pariser-Parr-Pople- Peierls model is a realistic and accurate model of 7r-conjugated 
systems, which includes the key features of long range electron-electron interactions and 
electron-lattice coupling. The Hamiltonian for an N site chain with open boundary condi- 
tions is defined as 

N-l i jV-1 N-l 

+ UJ2 (n iT - I) ( nil - + £ Vyim - - 1), (l) 

i=l v z/ v z/ <ij> 
where, < ij > indicates all pairs of sites, tg = (t + ^) and 

T e =W( c l+i*ci« + h.c.) (2) 

Z a 

is the bond order operator of the £th bond. We use the Ohno function for the Coulomb 
interaction: 

= U/y/l + fh*,, (3) 

where (3 = (f//14.397) 2 and bond lengths are in A. The dimensionless electron-phonon 
coupling constant, A, is defined by 

where K is the elastic spring constant (estimated to be 46 eV A~ 2 from Raman analysis of 
C-C stretching modes in trans-(CH) 2 ) p5| , and a relates the actual distortion of the £th. 
bond from equilibrium, Srg, to A^: 

5r e = A e /2a. (5) 

We take the undistorted chain to lie along the x— axis, with the bonds oriented at 30° to 
this axis. Then, for fixed bond angles, the distorted chain coordinates are defined as: 



4a 



ito = vS-iEA*(-l) m , (6) 



4 « l=i 
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where 

o \/3 , . 
x ij = ~y^\3 

and 

,o . 



Vij= °> if \j ~ *l even ( 7 ) 
= otherwise. 



ao (= 1.40A) is the undistorted C-C bond length. 
The force per bond, is 



f - 9{H) (R) 



Using the Hellmann-Feynman theorem this can be re-written as, 



f e =-2a[-^- + T-(T e )) (9) 



/ 



- S 2a(l + /?r&) 3 / 2 + ^"<^ - " 1}) J • 

The prime over the sum indicates that the sum runs over all pairs of sites which span the £th. 
bond. The contribution to the bond force from the Coulomb interaction is small compared to 
the kinetic energy term: the value of the Coulomb force from the nearest neighbor density- 
density correlator is approximately one tenth of the kinetic term. Moreover, the density- 
density correlator alternates in sign and drops to less than one tenth of the nearest neighbor 
density-density correlator, so the sum over all bonds is also small. Table II shows the 
correlator for up to five nearest neighbors. We therefore only include the nearest neighbor 
density-density correlator in the evaluation of the distorted geometry. (However, the full 
distorted geometry is used in the evaluation of the Coulomb interaction, Eqn. (3).) 

Using this approximation, and setting fe = 0, the self-consistent equation for the equi- 
librium Ai is: 

A f =y^)(TO-r-c> ), (io) 

where, 

c ' = iM^Mr <( "'" 1)(% '" 1)) ' 

We observe that, since the nearest neighbor density-density correlator is negative, the 
Coulomb interactions tend to increase the bond dimerisation. 

The calculations were performed for fixed chain lengths, which is enforced by setting, 



Y N-l 



N 



T £«7»-C*ao). (12) 
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To complete our discussion of the model we turn to its parametrisation. There are three 
parameters in the model: to, U and A. An optimal parametrisation for i and U was found in 
|T5| by fitting the Pariser-Parr-Pople model to the excited states of benzene. Assuming that 
this parametrisation is transferable between all 7r-conjugated systems, we use them here, 
and set to = 2.539 eV and U = 10.06 eV. The remaining parameter, A, is found by fitting 
the vertical energies of the i x B~ and 2 1 A£ states, calculated from the Parsier-Parr-Pople- 
Peierls model, to the 6-site linear polyene [IE |. This gives A = 0.115. Finally, using K = 46 
eV A" 2 implies a = 4.593 eV A -1 . 



III. SOLUTION OF THE PEIERLS MODEL 



As originally recognised by Pople and Walmsley the low lying excitations of the 
dimerised even N site chain correspond to the creation of two mid-gap states. These exci- 
tations are associated with localised geometrical structures which lead to a reversal of the 
lattice dimerisation, and were subsequently termed solitons. The defect states repel from 
each and are repelled from the ends of the chain. Thus, they reside at approximately N/4 
and 3N/4 along the chain (as may be seen in Fig. 5(b)). Fig. 1 shows a schematic energy 
diagram of the molecular orbitals and defect states, while Fig. 2 shows the energies of the 
l 1 ^; and 2 1 A+ states function of inverse chain length. It is clear that the first excited 
even parity state lies above the odd parity state. However, in the long chain, continuum 
limit, these states are degenerate, with energy 4A /7r = 0.12 eV, using A = 0.115, to = 2.539 



eV and 27 



8t exp 



(13) 



This gap is only a fraction of the experimentally measured gap of approximately 2.0 eV [28 
While a larger optical gap can be obtained by increasing A and to, the energetic ordering of 
the low lying states would still be incorrect. As we see in the next section, it is electronic 
interactions which primarily open the optical gap, and reverse the energetic ordering of the 
states. Furthermore, electronic interactions significantly modify the soliton structures, as 
we show in section VII. 

The non-interacting limit enables us to make a comparison between the DMRG methods 
and the exact calculation. In Fig. |3] the energy difference between the exact results and 
DMRG calculations is shown for the l l B~ and 2 1 A+ states. We see that for both states 
the accuracies of the infinite and finite lattice algorithm calculations are close, so that both 
methods can be used in the actual calculations. The accuracy is better for the l l B~ state, 
but even for the 2 1 A+ the error is about 0.002 eV for the 50 site chain in the infinite 
lattice algorithm calculation. Other DMRG convergence tests, confirming the validity of the 
method, were presented in |T6 . 



IV. SOLUTION OF THE PARISER-PARR-POPLE MODEL 

The uniform chain in the limit of only on-site Coulomb interactions is described by 
the Hubbard model. At half-filling, the spin excitations are gapless in the infinite chain 
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limit, whereas the charge excitations are gapped. Even though the Pariser- Parr- Pop le model 
contains long range interactions, the spin excitations still appear to be gapless in the uniform 
chain, as shown in Fig. 2. The 2*^4+ state is also gapless, confirming the interpretation of 
it as a pair of bound magnons. The optical gap (E(l 1 5~)) extrapolates to approximately 
1.6 eV, and is excitonic, lying approximately 1.0 eV below the charge gap for long chains. 
As discussed in section I, the energies of the 2 1 A+ and l x B~ states in polyacetylene thin 
films are believed to be at approximately 1.0 — 1.8, and 2.0 eV, respectively. Approximately 
0.3 eV should be deducted from the calculated l 1 ^" energy to account for solvation effects 
|29| |, indicating that the undimerised Pariser-Parr-Pople model underestimates the optical 
gap by approximately 0.7 eV and the 2 1 A+ energy by up to 1.8 eV. 



V. SOLUTION OF THE PARISER-PARR-POPLE-PEIERLS MODEL 



Sections III and IV indicate that neither electron-lattice coupling nor electron-electron 
interactions alone are sufficient to explain the low energy excitations of polyene oligomers. A 
pure electron-phonon model predicts degenerate l l B~ and 1 3 B+ states with the 2 X A+ state 
lying above them, while a pure electron interaction model underestimates the optical gap, has 
gapless spin excitations and does not lead to a dimerised chain. We now turn to the DMRG 
solution of the Pariser-Parr-Pople-Peierls model. We note that an infinitesimally small 
electron-phonon coupling will open a gap in the spin excitation spectrum for all electronic 
interaction strengths. 

We first calculate the ground state energy and lattice geometry. The normalised stag- 
gered bond dimerisation is defined as, 

^(-l)'(^, (14) 

where i is the average value of ti in the middle of the chain. 5 = 0.102 in the center of the 
chain. Using a = 4.593 eV A -1 , this implies that the bond length alternation of the ground 
state in the middle of the chain is 0.056 A, in close agreement with the experimental result 
of 0.052 A |0|. 



Using the ground state geometry, the vertical energies (that is, the energies of these states 
with the ground state geometry) (E v ) of the 1 3 B+, 1 X B~ and 2 X A+ states are calculated. 
These, as well as the relaxed energies (E°~°), are shown in Fig. 4(a) as a function of inverse 
chain length. The vertical energy of the 2 X t4+ state lies approximately 0.3 eV above that 
of the l l B~ state in the long chain limit The relaxation energy of the ^B~ state is 



modest, being approximately 0.2 eV for 102 sites. By contrast, the relaxation energies of the 
1 3 £>+ and 2M+ states are substantial, being approximately 0.8 eV and 1.5 eV, respectively, 
and converge rapidly with increasing chain length. The energy of the relaxed 2 X A+ state 
lies 1 eV below that of the l 1 -^ state. We see in section VII that this strong relaxation is 
associated with a large distortion of the ground state structure. 
In Fig. 4(b) we plot the charge gap, 

E(N+1)+E(N -1) -2E(N), (15) 

and the energy of the ^B~ state. In the long chain limit the charge gap represents the 
energy of an uncorrelated electron-hole pair, and therefore represents the band edge. The 
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relaxation energy of the charge gap is roughly double that of the l 1 -^ state. This is to be 
expected, as the two charges form independent polarons, whereas the excitonic 1 X B~ state 
forms a single polaron, as discussed in section VII. We see that the single chain binding 
energy is 2.4 eV. However, the unbound pair is strongly solvated (ca. 1.5 eV), while the 
exciton is more weakly solvated (ca. 0.3 eV) [P5fl . This implies that the bulk binding energy 
of the 1 X B~ state is ca. 1 eV. 

The experimental values of i? 0_0 (l 1 i?~) and E°~°(2 1 A'^) for short polyenes are also shown 
3~2fl . The 2 1 A+ values are in excellent agreement with our calculation. The 1 X B~ values 



are approximately 0.3 eV lower than our predictions, which is approximately the reduction 
expected by the solvation of the chains in solution Thus, for short polyene oligomers, 
the optimised parametrisation of the Pariser-Parr-Pople-Peierls model gives remarkably good 
results. 

Kohler has analysed the experimental results for N = 6 — 16 [32|. For the 2 1 A+ state 
the empirical relation, 

E°-°(2 1 A+) = 0.96 + 20.72/iV, (16) 



was derived. This relation appears to confirm the work of who find a 2 1 A+ state at 
1.1 eV in thin films. However, there is no particular reason why a linear extrapolation in 
1/N is valid. Our calculation for the Pariser-Parr-Pople-Peierls model shows a significant 
flattening off of the 2 1 A+ energy for chain lengths of roughly 30 or more sites. The calculated 



converged energy of 1.74 eV is in agreement with |[13| . 

This rapid convergence of energy with chain length is in contrast to both the Pariser- 
Parr-Pople and Peierls models. In the Peierls model the excitation energies are gapped, 
but the deviation from 1/N behavior is only evident for long chains (ca. 100 sites). In 
the Pariser-Parr-Pople model a deviation from 1/N behavior is only evident in the long 
chain limit for the 1 X B~ state and the charge gap. In the Pariser-Parr-Pople-Peierls model, 
however, states which form pronounced solitonic structures, such as the 2 1 A+ and triplet 
states (as discussed in section VII) self-trap once the chain length exceeds the size of their 
solitonic structures. It is possible that this self-trapping is a consequence of the adiabatic 
treatment of the lattice, and that a full treatment involving quantum phonons would change 
this prediction. 

Our understanding of self-trapping - and its validity or otherwise - is complicated by the 
discussion of the 1 X B~ state energy. Again, an empirical relation, 

^(l 1 ^-) = 2.01 + 15.60/iV, (17) 

was derived by Kohler, which is in good agreement with the thin film result. Our calculated 
value of 2.74 eV is too high, even when solvation effects (ca. 0.3 eV) are deducted. Once 
again, the 1 X B~ state is self-trapped, and the possible relaxation by lattice fluctuations 
would lead to a better agreement. 

However, since the phonon frequency of ca. 0.2 eV is so small compared to the electronic 
energy scales, any corrections to the adiabatic limit are expected to be small, so we need to 
consider other possible reasons for the discrepencies in the long chain limit. One source is 
the possible renomalisation of the Pariser-Parr-Pople-Peierls model parameters in the long 
chain limit; another is a-electron screening. 
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VI. PHOTO-INDUCED ABSORPTION 



The photo-induced absorption spectrum of a system, obtained while it is being pumped 
at an energy above the optical gap, gives an insight into the excited states of that system 
| 12|1 . Typically the system is pumped at 2.4 eV, and photo-induced absorption peaks are 
observed at 0.43 eV and 1.35 eV. 

The higher energy peak is believed to intrinsic and has been ascribed to a bound soliton- 
anti-soliton pair ]33|]. A possible interpretation is that excitations to states above the vertical 
l 1 ^?" state decay non-radiatively to the 2 1 A+ state, which subsequently relaxes. The photo- 
induced absorption is then a vertical transition from the relaxed 2 X A+ state to a l B~ state. 
We find that the energy of the l 1 -^ state in the relaxed geometry of the 2 1 A+ state lies 
1.3 eV above the 2 1 A+ state for 102 sites. However, the dipole moment is weak, being 
only 0.16(/x) 1 i_ B - (where ([i)ii B - is the dipole moment between the ground state and 1 X B~ 
state). A second possibility is that it is a triplet-triplet (T — ► T*) transition. We calculate 
this transition energy to be 2.8 eV, while the dipole moment is 0.96(/i) 1 i B - . Since the excited 
triplet (T*) is a triplet-exciton (as opposed to a spin-density wave excitation) at high energy, 
it is reasonably to assume that it will be strongly solvated, reducing this transition energy 
by as much as 1 eV. Thus, a triplet to triplet transition is a possible explanation for this 
absorption. 

We calculate the transition energy between the lowest polaron state and the first dipole 
connected excitation to be 0.45 eV at 102 sites, and the dipole moment is 0.88(/i) 1 i B -, 
suggesting that this is the origin of the lower peak. 



VII. SOLITON STRUCTURES 

In Fig. 5(a) we plot, as a function of bond index from the center of the chain, the 
normalised staggered bond dimerisation, 5g, Eqn. (14). We note that the 1 3 B+ and 2 1 A+ 
states undergo considerable bond distortion, whereas the l l B~ state shows a weak polaronic 
distortion of the lattice, similar to the distortion associated with a doped charge. In [16[] we 
showed that the 1 3 5+ and \ X B~ states fit a 2-soliton form 0, [Q], ]35J, whereas the 2 1 A+ 
state fits a 4-soliton form. The bond distortions of the non-interacting limit (the Peierls 
model) are plotted in Fig. 5(b). A comparison between these plots illustrates the role played 
by the electronic interactions in modifying the non-interacting picture: 

1. The dimerisation in the ground state is enhanced by a five-fold factor, in qualitative 



agreement with 19] and 20 



2. The 1 1 B U state evolves to an exciton-polaron, in agreement with [|Hj 



3. The 2M+ state, owing to its strong triplet-triplet contribution, evolves to a 4-soliton 
solution, in agreement with ||. 

Further insight into the electronic structure of polyenes and its relation to their geometry 
can be obtained from the spin-spin correlation function, defined as, 

S l = -{SfS z N+1 _ l ). (18) 
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This function measures anti-ferromagnetic correlations between sites symmetrically situated 
with respect to the center of the chain. As the correlation function shows unimportant 
oscillations between even and odd site indices i, we use the symmetrized function: 

Sj = 2^ S h( N ~i) + ^(iv-^+i)' ( 19 ) 

j — 0, 4, 8, N — 2, which measures the correlations between pairs of doubly-bonded sites, 
with j being the distance between them. 

The spin-spin correlation functions, calculated in the ground state geometry, are shown 
in Fig. 6(a). They show a monotonic decay for the correlations in the l 1 ^ and ^B' states, 
but in the 2 1 A+ state there is a small minimum at j = 8 and a maximum at j = 16. This 
behavior of the spin-spin correlations in the 2 l A+ state becomes clearer when we calculate 
it in the relaxed geometry for this state. Here, the correlation function of the 2 1 A+ state, 
shown in Fig. 6(b), has a strong minimum at j = 8, where it changes sign, and a maximum 
at j = 20. These features strongly confirm the triplet-triplet character of this state. By 
comparing Fig. 6(b) to the soliton structure shown in Fig. 5(a), we see that the unpaired 
spins correspond to the positions of the geometrical solitons. 

VIII. DISCUSSION AND CONCLUSIONS 

We began this investigation of the electronic and geometrical structure of linear polyenes 
by performing seperate studies of the U — Peierls model and the A = Pariser-Parr-Pople 
model. These studies show that these two limits predict quite different low-lying excitations. 
The Peierls model predicts mid-gap states associated with geometrical defects. The dipole 
forbidden 2 1 A+ state lies above the degenerate singlet and triplet 1B U states. In contrast, the 
Pariser-Parr-Pople model predicts gapless (or very small gapped) triplet and 2 1 A+ states, 
with the 1 X B~ state lying above them. 

When these two models are combined in the Pariser-Parr-Pople-Peierls model we see the 
effect of the inter-play of electron-electron and electron-phonon interactions. The lowest 
lying triplet (l 3 B^) is a soliton-antisoliton pair; the lowest lying singlet (2 1 A^) is an even- 
parity pair of soliton-antisoliton pairs, owing to it being a bound pair of triplets; and the 
lowest optically allowed state (^B') is an exciton-polaron. The soliton positions in the 
2 1 A+ state is confirmed by the spin-spin correlation function. Electron-electron interactions 
play the dominant role in opening the optical gap and dimerising the lattice. 

We find that the relaxation energy of the l 3 B+ and 2 1 A+ states are substantial, whereas 
that of the 1 X B~ state is modest. The vertical energy of the 2 l A+ state lies above that of 
the l 1 -^" state, but the relaxed 2 l A+ state lies ca. 1.0 eV below that of the l l B~ state. 
The role of electron-electron interactions are crucial and subtle in determining these relative 
positions. A larger electron-electron interaction leads to a more dimerised ground state, and 
this tends to raise the vertical energy of the 2 l A+ state relative to that of the l l B~ state. 
However, a larger electron-electron interaction also leads to a larger relaxation of the 2 1 A+ 
state energy compared to that of the l l B~ state, leading to a reversal of their energies. 

For short polyenes we find good agreement with experimental values. However, in the 
long chain limit the results (at least for the 1 X B~ state) become more qualitative. The 
experimental uncertainty in the position of the 2 1 A^ state means that we cannot be sure of 
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the validity of our prediction. However, if we assume that ca. 1.0 eV is the correct relaxed 
energy of the 2 1 At state, then our predictions are between 0.5 to 1.0 eV too high. In section 
V we discussed some of the possible origins of these discrepencies. They include, the neglect 
of lattice fluctuations in the adiabatic treatment of the lattice, the possible renormalisation 
of the 7r-model parameters in the long chain limit, and the neglect of the a-bond screening. 
We would expect that as a molecule gets larger the 7r orbitals will become more extended, 
as they mix with other orbitals. This will reduce U and a (and hence A), and increase to, 
thus reducing the excitation energies. Work is currently in progress to study these affects. 
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TABLES 



TABLE I. The classification of the relevant states. 



State L4+ LE^ 3 B+ 

Spatial Inversion Symmetry + — — 

Spin 1 

Particle-hole symmetry + — + 

Character Covalent Ionic Covalent 



TABLE II. The density-density correlator as a function of distance 



3 (Jjh - Ife+j - !)) 

1 -0.308 

2 +0.002 

3 -0.021 

4 +0.004 

5 -0.011 
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FIGURES 



FIG. 1. Energy level diagram for the key low-lying states in the non-interacting limit. 

FIG. 2. Transition energies for the ^B~ (squares), 2 1 J 4+ (diamonds), 1 3 -B^ (triangles) states 
and charge gap (circles) as a function of inverse chain length for the U = Peierls model (dashed 
lines and open symbols) and the A = Pariser-Parr-Pople model (solid lines and filled symbols). 
(In the Peierls model the l 1 ^" and 1 3 -B^~ states are degenerate.) 

FIG. 3. The difference between the exact calculation of the 2 1 ^4+ (diamonds) and ^B~ 
(squares) states in the non- interacting limit, and the results of DMRG calculations in the infinite 
and finite algorithms. Solid lines correspond to infinite lattice algorithm results, dashed lines to 
the finite lattice algorithm. 

FIG. 4. (a) Transition energies for the l l B~ (squares), 2 1 ^j" (diamonds) and 1 3 -B+ (trian- 
gles) states as a function of inverse chain length. Vertical/relaxed transitions are indicated by 
dashed/solid lines and open/solid symbols. Experimental results for the relaxed ^B~ (x) and 
2 1 ^j" (+) state energies for polyenes in hydrocarbon solution |^2[ . (b) Transition energies for the 
^B~ state (squares) and charge gap (circles) as a function of inverse chain length. 

FIG. 5. (a) The geometries (normalised staggered bond distortion bi as a function of bond 
index t from the center of the lattice) of various states of the Pariser-Parr-Pople-Peierls model: 
(crosses), 1 1 B U (squares), l 3 B+ (triangles) 2 1 A+ (diamonds) and polaron (circles), for the 
102 site system, (b) The same as (a) for the U = Peierls model. 

FIG. 6. Spin-spin correlation functions for l 1 ^^ (solid squares), 2 1 A^" (solid diamonds) and 
l 1 !?" (empty squares) states, (a) In the relaxed l 1 ^^" geometry, (b) in the relaxed 2 1 A+ geometry. 
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